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A Quantum Energy Inequality (QEI) is derived for the massive Ising model, giving a state- 
independent lower bound on suitable averages of the energy density; the first QEI to be established 
for an interacting quantum field theory with nontrivial S-matrix. It is shown that the Ising model 
has one-particle states with locally negative energy densities, and that the energy density operator 
is not additive with respect to combination of one-particle states into multi-particle configurations. 
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In relativistic quantum field theory, the Hamiltonian 
is a positive operator in all inertial frames of reference; 
this is the content of the spectrum condition. By con- 
trast, it is impossible for nontrivial local averages of the 
energy density to be positive operators in any quantum 
field theory obeying standard assumptions [T]: typically, 
the expectation value of the energy density at any point is 
unbounded from below with respect to the state [2] . This 
state of affairs demonstrates a fundamental incompati- 
bility between quantum fields and the energy conditions 
usually assumed in classical general relativity, and which 
are the essential input for results such as the singular- 
ity theorems of Penrose and Hawking [3H5] , the positive 
mass theorems [BHS]) and Hawking's chronology protec- 
tion results [TU], among many others. 

Nonetheless, various models of quantum field theory 
obey local remnants of the spectrum condition called 
Quantum Energy Inequalities (QEIs), which provide 
lower bounds on expectation values of the energy density 
when averaged along a timelike curve or over a spacetime 
region. Their study originates from Ford's insight [11] 
that quantum field theory could produce observable de- 
viations from the second law of thermodynamics unless 
there were mechanisms to constrain negative energy den- 
sities (or fluxes) . QEI bounds place severe constraints on 
the extent to which quantum fields can support exotic 
spacetime geometries [T^ |T3] ; moreover, weakened clas- 
sical energy conditions inspired by QEIs can be used to 
prove singularity theorems |14j . 

QEIs have been established for free (minimally cou- 
pled) Klein-Gordon [15H2I], Dirac [22H2g, Maxwell |26l 
\Tf\ and Proca fields [37] in both flat and curved space- 
times, for the Rarita-Schwinger field in Minkowski 
space [25], and for the whole class of unitary positive- 
energy conformal field theories in two-dimensional 
Minkowski space [29] (generalising a special case [l6]). 
In all these cases, the energy density is bounded from 
below on the class of physically acceptable states, if it 
is smeared against a positive test function over a re- 
gion or curve of nonzero temporal extent. Even within 
the setting of free fields, however, the nonminimally cou- 
pled Klein-Gordon field provides an example where only 
a weaker type of QEI holds [30" - the lower bound is 



no longer state-independent, but exhibits dependence 
on the energy scale - and it is expected that this be- 
haviour would be typical for interacting quantum field 
theories [3T] . Analogues of these energy-dependent QEIs 
exist in general quantum field theories for observables 
arising in operator product expansions of 'classically pos- 
itive' expressions [32]; however, a direct connection to the 
energy density is lacking. 

In this paper, we derive a QEI for the continuum Ising 
model |33| . a quantum field theory of interacting mas- 
sive spin-0 bosons in two-dimensional Minkowski space. 
The Ising model is the simplest integrable quantum field 
theory (other than free models) with a factorizing scat- 
tering matrix - the two-particle 5'- matrix is 5*2 = —1. It 
is closely related to other integrable models such as the 
sinh-Gordon theory. Our result is the first QEI to be ob- 
tained for a theory with a nontrivial S'-matrix. Although 
the scattering theory is rather simple, the energy density 
of the Ising model differs in essential ways from that of 
the free scalar field. As we will show, there are single- 
particle states of the Ising model that have locally nega- 
tive energy density. In a free theory one could combine 
such single-particle states to obtain multi-particle config- 
urations with arbitrarily large negative energy densities 
(cf. [30]), thus excluding the existence of a state inde- 
pendent QEI. However, the energy density of the Ising 
model is not additive with respect to tensor products of 
single-particle states; a signature of the interacting na- 
ture of the theory. It is intriguing that the interaction is 
responsible for maintaining the QEI in this sense. 

Although it is a theory of bosons, the Ising model is 
conveniently defined on a fermionic Fock space. Our ar- 
gument will make essential use of this underlying struc- 
ture, and indeed is based on arguments originally devel- 
oped in the context of the free Dirac field [3^24j . The 
QEI bound is given as the (— ) case of ([8| below; the 
{+) case applies to the free scalar field and is given for 
comparison. 

To start, we fix our conventions for the quantum 
Ising model of mass /i > 0, in a notation that makes 
parallels to the two-dimensional scalar free field clear. 
We work on the single particle space /C = L'^{R,d9), 
where the rapidity 6 is related to two-momentum p by 



p{9) — /i(cosh0,sinh6'). The two models (with (+) for 
the free field and (— ) for the Ising model) are defined on 
the symmetric, respectively antisymmetric, Fock space 
'H± over /C, with vacuum vector denoted as fl. On 'H±, 
we have the usual action of creators and annihilators 
a±{9), a±{r]), which fulfill canonical (anti-)commutation 
relations, 

a±iv)o-iiO) T aiie)a±{r^) - 5{0 - rj)l. 

Spacetime symmetries, i.e., translations x ~ {t, x), boosts 
A, and the space-time refiection j, act on 'H± hy 

C/±(x,A)al(0i)---4(0„)^ = 

g*(p(ei)+...+p(9„)) V^(0^ + A) • • • aliOn + X)n, (1) 



fermionic Fock space H- by 



u±{])al{ei) ■ ■ ■ alie„)n = alien) 



J 



J 



h)n, 



(2) 



with U±{j), the PCT operator, extended antilinearly. 

We now describe the basic observables of the model. 
Following Lechner |34j (though with slightly different 
conventions) , we define quantum fields 4>± , 0^ as 



1 



de fe'PW"<a^(6') +e-'PW-^a±(6l)V 



'/'±(x) 



'/''±(x)-C/±(j>±(-x)(7±(i). 

These are covariant under the symmetry operations M, 
^ , with space-time reflection exchanging (/)± and (f>'±. . Of 
course, (/>+ — (f)\ is just the usual local free scalar field. 
But 0_ is an interacting and nonlocal field; rather than 
local commutation relations, we have [0_(x), 0'_(y)] — 
if (x — y)^ < and x^ < y^ [we say x is to the left of y]. 
An operator A is considered to be localized at x (or in a 
region O) if 



i_(y),A] = O=[A0-(z)] 



(3) 



whenever y is to the left of x (or O) and z is to the 
right. By abstract arguments, a large class of such local 
observables exists 34]. 

Let us recall the scattering theory of our models. The 
case 0+ is trivial of course. As mentioned, (j)- defines the 
(interacting) massive Ising model with nontrivial scatter- 
ing matrix |34| . Both the incoming and outgoing states 
are bosonic Fock states, i.e., we can identify Hin and Hout 
with H+. The incoming M0ller operator is then given by 
Vin : H- —7- H+j 



Vi^aUoi) ■ ■■aUOn)n = (l[e{e,-0,))alie,) ■ --alie. 



l<] 



and the outgoing M0ller operator V^ut by a similar for- 
mula, but with the argument of each e negated. The 
S'-matrix is S = Mn^^out = (-l)^+(^+-i)/2, where N+ is 
the bosonic number operator on "H^. 

A crucial observation, for our purposes, is that a free 
Majorana field ip = ('01:^2)'^ can be defined on the 



■01,2 (x) 



^Q g..(-2±l)/4g^p(e).x±e/2^t_ (^) ^ ^^ 



(here -I- for the case 1 and — for the case 2). This Ma- 
jorana field is covariant under f/_(x. A) as defined above, 
but not under U-{j); its associated PCT operator is fun- 
damentally different. By analogy with '35, Sec. 6], 
fulfills 



{0_(y),V(x)} = O=[0(y),0'_(x)] 



(4) 



if y is to the left of x. That is, V'(x) is not a local observ- 
able of the interacting theory (cf. (Is])). However, from 
(HI), all even polynomials in 0i 2 have this property. This 
applies in particular to the energy density 

TOO(x) = '--.^{xfdtHx) - {dt^{x)fi^{x): . 

Since the energy-momentum operators of the Majorana 
field coincide with those on "H- by (fTl), T^^ is also the 
energy density of the interacting Ising model; nonethe- 
less, we emphasize that the Ising model is distinct from 
the free Majorana theory. 

For completeness, we recall that the energy density of 
the free Bosonic scalar field 0+ on H+ is 

Tl'{x) - lidtMx))' + {d.Mx)? + /i'0+(x)': . 

After some computation, the energy densities of both the 
free and interacting Bose models take the form 



Tf 



(x) ^^J dddr](^F±ie,ri,x)aii9)ai{r,) 

+ 2F±{9,r, + 7Tt,x)ali9)a±{7j) (5) 

+ F±{d + ni,!] + TTi,x)a± {9)a± (77) 



with 



F+iCx) = -^sinh^ ^L±^e*(p(Ci)+P(C.))-x^ 
27r 2 

F_(C,x) =^zsinh '^^~^^ F+(C,x). 



(See [3S1 [37] for the general theory of these expan- 
sions.) Using methods from [371 Sec. 9.1], one can 
see from the analyticity structure of F± that T™(/i) = 
\Q J dth{t)Tj°{t,x) is a closable operator, and a local ob- 
servable, for any a; € K and Schwartz test function h. 

In the Bosonic free field situation, it is well known (see 
e.g., [35]) that T™ has nonnegative expectation value in 
all single particle states, or more generally, states with 
sharp particle number, while more general superpositions 
(such as that of the vacuum and a two-particle state) 
can yield negative expectation values (see, e.g., [T]). By 
contrast, we will now exhibit single-particle states in the 



Ising model with negative energy density at the origin, 
reminiscent of the situation for free Dirac fields [38] . 

These states are essentially superpositions of two plane 
waves. More specifically, let us choose a nonnegative real- 
valued Schwartz function h with J h{9)dd — 1. We set 
(f> :— a_{ipa,p,'y)^ = J dO ipa,p,'y{9)aL{9)Cl, where 

with parameters a > 0, /3, 7 G M and normalization con- 
stant c„,;3^^ > 0. We will show that ($,r™(0)$) < 
for a suitable choice of the parameters. To that end, we 
compute from ([5|, 

($,rO"(0)a>) = '^Clp,4lc + Ja,7/3 + ^a,^/?'), (6) 

where we denoted 
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/a = / d6ldr; /la (6l)/ia (77) cosh^ ^ cosh ^, 

X„^^ = f d6dr]haie)hair]) cosh^^^^^^ cosh ^. 

The right hand side of (|6| is negative for some /3 if the 
polynomial /„ -|- Ja^-fP + Ka^^P'^ has two real zeros, that 
is, if J'^_^ > AlaKa^f This inequality holds for small a 
if it holds in the limit a — )■ 0. Noting that ha{9) -^ S{0) 
in this limit, we obtain 



Figure 1. (Color online) Energy density per particle in the 
Ising model, computed in the multi particle states $„_ along 
the time axis, with parameters a = 0.5, /3 = —0.04, 7 = 5, 
h{e) = exp(-6|2). 



ri($i+,r™(x)<i>i+). However, the same relation does not 
hold in the interacting situation: In general. 



$„_,Tl!"(x)<&„_> ^ n($i_,r!'"(x)$i_ 



(7) 



Iq — ^, Jo.- 



2 cosh' 



3 7 



Kqj = cosh 7. 



We can deduce this from our other results: We saw above 
that ($i_,T™(x)$i_) is negative in some examples, and 
therefore at large n, equality in ([7| would be in contra- 
diction to the QEI (l8| that we will establish below. 

But let us give a direct argument for (l7| . It is useful to 
note that Vinal{e)V*^ = a\{e)M{e), where M{e) is the 
multiplication operator 

M(0)a^77i)---a;(77„)fi = 



j=i 



Using this relation and its adjoint, it is straightforward 
to compute from ([5| that 



2 



The condition Jq ,^ > AIqKq^^ then becomes 
(1 -I- cosh 7)^ > 8cosh^7, 

which is in fact fulfilled for sufficiently large 7. With 
these choices, we achieve ($, T™(0)$) < 0. For an exam- 
ple see Fig. [1] (solid line), where ($,r™(i,0)$) is plotted 
for a suitable choice of the parameters. 

Given any normalized single-particle state wave func- 
tion if £ K., we may form multi-particle states by taking 
tensor products in the 'in' Hilbert space and applying 
the inverse M0ller operator. For the free model, this 
yields n-particle vectors $„+ := {nl)~^^^a\^{ip)"D, e 1-1+, 
while the corresponding states in the Ising model are 

*„_ :=yi;$„+eH_. 

Now, the total energy, given by the Hamiltonian H± 
(the generator of time translations), is additive in the 
sense that ($„±,i?±$,i±) = n{^i±,H±^i±). (For H+ 
this is evident from (fTj), and for iJ_ it follows since 
VinH^V^^ = H+.) Furthermore, the energy density in the time axis for the same choice oi if = fa,^,^ as above, 
the free model is also additive: ($„_|_,T5'°(x)$„+) — This behaviour refiects the interaction in the Ising model. 



($„_, rO"(x)$„_) = n ^ / dOdricosh' ^^ cosh 
2n J 

X ^^{tj)L^{0, ,;)"-ie»(pW-pW)"< 
where 



L^{e,Tj):= d\\^(\)\ 



X)e{r,-X). 



The factor L^, which does not occur in the free case, pre- 
vents additivity of the energy density. This nonlinearity 
is also apparent in Fig. [l] where the energy density per 
particle, p{t) := i($„_, r™(i, 0)$„_), is plotted along 



Finally, we turn to the derivation of the QEI bounds on 
the energy density. In fact, the free scalar field (+) and 
the Ising model (— ) obey closely related QEIs, namely 



dig(i)^(<l>,rr(i,x)$)> 



1 

'4^ 



diouj^\giu;)\^Q±{-), 

(8) 
for any real- valued smooth compactly supported function 
g and all sufficiently regular normalized states $, where 
the dimensionless functions Q± : [l,oo) — > M+ are 



Q±{u) = Vl - u-^ ± u-^ log(u + Vu^ - 1) (9) 

and obey Q±(l) = 0, \imu^+oDQ±{u) = 1, while the 
Fourier transform is g{uj) = J dt g{t)e'^'^* . 

The QEI for the free Bose field (+) has been known 
for some time [15] and holds rigorously for all Hadamard 
states $ [5T] ; our goal is to establish the QEI for the Ising 
model (— ). As we have already seen that its energy den- 
sity coincides with that of the free Majorana field, the 
QEI is exactly the same as for the latter theory. This 
has not been computed before, but the argument is suf- 
ficiently similar to treatments in [531 IM] that we only 
sketch it here. For i^ e M, let R^ and 5*^ be continuous 
one-parameter families in /C = i~(M, d^), and define for 
each V the (bounded) operator 



de (R^{0)a^{d) + SAe)al{e) 



O, 



on T-L-. Elementary use of the CARs shows that 

OlO, - \\S,\\H = -O^Ol + \\R,\\H =: X,, 

where || • || is the norm in K.. As uO]^0^ (resp., —vOyO'l) 
is positive scmidcfinitc for v >{) (resp., u <Q),we have 



dv 



v(^.X,M > 



dv 



^{\\R-A. 



\SA 



(10) 
for all normalized quantum states $. For our application, 

— qiv — n cosh 6) cosh - , 

n 

q{v -\r fJ. cosh 9) sinh - , 
47r' 2 

for any real- valued g e C^ (M) . Using the identity 

' di' 




{uj + J)g'^{J -cj) = - 



-vg{v + uj)g{v + cj') 



(a mild rewriting of Eq. (2.17) in [23]), the left-hand side 
of ([lO|) becomes, after a computation. 



dv 



v(^,X,M 



dtg{t)^{^,Tl\t,Q)^). 



It is expected that this is rigorously valid at least for all 
<& € H- that are Hadamard states of the Majorana field 



(cf. the treatment of the Dirac equation in four dimen- 
sional curved spacetimes [25]). We also compute 



RHS of (101 



'47r2 



47r2 



OO /'OO 



du u I dO cosh 6\g{^co&]i 6 -\- h')Y 



(11) 



duuj \g{u)\ Q^{uj/n), 



where Q- was defined in (|9|. In more detail, the sec- 
ond equality in (111 is obtained by using evenness of 
the integrand to alter the inner integration region to 
[0, oo), then changing variables from {ly, 9) to (w, 9) where 
uj — 1/ + ^ cosh 9 (with a consequent change of integration 
region) and finally evaluating the 6 integral. For x — 0, 
the energy density of the Majorana and Ising models thus 
satisfies the QEI given as the (— ) case of (|8|, for all real- 
valued test functions g and for a large domain of $ G H- ; 
the result for general x follows by translation invariance. 
This bound is precisely half of the bound for the free mas- 
sive Dirac field in two dimensions obtained in ;24] using 
similar arguments, as might be expected. 

It is also worth considering the limit of these QEIs as 
/i — >■ 0, for fixed g, corresponding to the short-distance 
scaling limit of the theory [35]. In both cases. 



RHS of 



1 

47r 



dtwm' 



where we have used the fact that |5(u;)| is even, and the 
Plancherel theorem. But both the massless free scalar 
field and the massless Majorana field are conformal field 
theories, and obey sharp QEIs [2S1 Eq. (4.25)] 



digit)'' 



($,rf (t,a;)$) > 






dt\g'{t)\' 



for suitable normalized $, where C± arc the central 
charges of the left- and right- moving components: C+ = 
1 (free scalar field) and C_ — | (free Majorana), so the 
sharp bound is therefore tighter by a factor of 3/2, re- 
spectively 3, in these two cases. Accordingly, we do not 
expect our QEI to be sharp for fj, > 0. 

In summary, for the first time, a QEI has been derived 
for a quantum field theory with nontrivial S'-matrix. We 
hope this will prove to be a foundation for similar results 
on other integrable theories in two spacetime dimensions. 
Intriguingly, the interaction conspires to yield a QEI with 
a state-independent lower bound; it will be interesting to 
see whether this persists in other models. 

The authors thank Atsushi Higuchi for his perceptive 
questions and comments. 
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